Status of the differential transformation 

method 



C. Bervillier 

Laboratoire de Mathematiques et Physique Theorique, 
UMR 7350 (CNRS), 
Federation Denis Poisson, 
Q-| Universite Francois Rabelais, 

Pare de Grandmont, 37200 Tours, France 



5-H 



Of 



Abstract 

Further to a recent controversy on whether the differential transformation method 
(DTM) for solving a differential equation is purely and solely the traditional Taylor 
series method, it is emphasized that the DTM is currently used, often only, as a 
technique for (analytically) calculating the power series of the solution (in terms of 
the initial value parameters). Sometimes, a piecewise analytic continuation process 
is implemented either in a numerical routine (e.g., within a shooting method) or in 
a semi-analytical procedure (e.g., to solve a boundary value problem). Emphasized 
also is the fact that, at the time of its invention, the currently-used basic ingredi- 
ents of the DTM (that transform a differential equation into a difference equation 
of same order that is iteratively solvable) were already known for a long time by the 
"traditional" -Taylor- method users (notably in the elaboration of software packages 
-numerical routines- for automatically solving ordinary differential equations). At 
now, the defenders of the DTM still ignore the, though much better developed, stud- 
ies of the "traditional" -Taylor- method users who, in turn, seem to ignore similarly 
the existence of the DTM. The DTM has been given an apparent strong formal- 
ization (set on the same footing as the Fourier, Laplace or Mellin transformations). 
Though often used trivially, it is easily attainable and easily adaptable to different 
kinds of differentiation procedures. That has made it very attractive. Hence appli- 
cations to various problems of the Taylor method, and more generally of the power 
series method (including noninteger powers) has been sketched. It seems that its 
potential has not been exploited as it could be. After a discussion on the reasons of 
the "misunderstandings" which have caused the controversy, the preceding topics 
are concretely illustrated. It is concluded that, for the sake of clarity, when the DTM 
is applied to ODEs, it should be mentioned that the DTM exactly coincides with 
the traditional Taylor method, contrary to what is currently written. 
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1 Introduction 



The differential transformation method (DTM) of Pukhov [JJ [5] and Zhou 
[6] is frequently presented as a (relatively) new method for solving differential 
equation^]. Though based on Taylor series, it would be different from the 
traditional Taylor (or power) series method presented in usual textbooks as 
e.g., [7J. This distinction has been the object of a recent dispute [8f9] . Indepen- 
dently of whether this distinction was present in the original ideaJ 2 ! it is at 
least seemingly clearly expressed and often repeated since the second half of 
the 90's when the DTM has more systematically been used "to solve differen- 
tial equations" [I0]-[16]. F° r example, in [H] the "differential transformation 
technique", is presented "as an extended Taylor series method" and in [To] , 
where the differential equation referred to is 



one can read (in addition): 

"Tie differential transformation technique is one of the numerical methods 
for ordinary differential equations. It uses the form of polynomials as the 
approximation to exact solutions which are sufficiently differentiable. This is 
in contrast to the traditional high-order Taylor series method, which requires 
the computation of the necessary derivatives of f(x, t) and is computationally 
intensive as the order becomes large. Instead, the differential transformation 
technique provides an iterative procedure to get the high-order Taylor series. 
Therefore, it can be applied to the case with high order." [T5J p. 25] (see also 



That statement is important since then the idea that the "traditional" Taylor 
series method would require the explicit calculation of high-order derivatives of 
/ (x, t) is repeated in many articles on the DTM as if it would be a convincing 
reason to make a definitive distinction between the two methods. 



Email address: claude.bervillier@lmpt.univ-tours.fr (C. Bervillier). 

1 An extensive presentation of the DTM is given in section 12.21 

2 The most frequently cited original works on DTM are inaccessible to me because 
they are written either in Russian |4|5] or in Chinese [6]. Except those two refer- 
ences, I have systematically solely cited articles written in English and, hopefully, 
accessible to all. 
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Before going any further, it is worth recalling what is the Taylor series methocjjO: 



Definition 1 ((Formal or raw) Taylor series method) The Taylor series 
method consists in expressing the solution of (TJP as a power series expansion 
about the initial time t : 
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in which the derivatives 

at 



are such that is satisfied order by order 

to 

in powers of (t — to). As consequence, the Taylor coefficients of the expansion 
are completely determined once the initial parameter «o is fixed. 

Notice that this definition specifies neither the way the derivatives ^ ^ are 
calculated nor the convergent property of the series. Thus it is surely clear to 
all that the formal Taylor series method (described in textbooks such as [7]) 
requires only that the series of x (t) be obtained by any means whatsoever. 
Consequently, it is obvious that if the DTM is a "procedure to get the high- 
order Taylor series" pT5lfT6] . then the method used to solve the differential 
equation is not the DTM but the Taylor series method (at least formally). 
Consequently to really appreciate the above quotation of [15], it is necessary 
to understand what actually is the "traditional" Taylor series method. No 
doubt that for the DTM users, the naive explicit computation of high order 
derivatives of f{x,t) is an integral part of the "traditional" method. One of 
the object of the present article is to show that this allegation is false. 

In fact, as one may anticipate from the presence of the terms "technique" and 
"numerical method", the Taylor series method mentioned in the above quo- 
tation of [15], though not explicitly defined, is not the formal ond^l It would 
be possible that it is a problem of semantics which has generated misunder- 
standing. 

Those words (technique, numerical) have wittingly been used to qualify the 
DTM. Indeed, some authors would like to see the DTM not as a formal 
method, but exclusively as a new numerical approach: 

"... the construction of power-series solutions has been generally thought of 
as an analytic tool and not as the basis for numerical algorithms. This is 
changing, and algorithmically constructing power-series solutions to ordinary 



3 For the sake of a convenient writing, extensive descriptions of the "traditional 
Taylor series method and of the DTM are postponed to section [2j 

4 Which is the best known formulation of the Taylor method. 
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differential equations is gaining in popularity. This is now often called the 
Differential Transform(ation) Method (DTM)." pi] 

However, if this view corresponds well to Pukhov's initial aim of "consider [ing] 
the problem of feasibility of constructing computer-specialized procedures ori- 
ented toward automatic solution of Taylor equations" [T] , the actual use of 
power-series solution as the basis for numerical algorithms is not at all new, 
see, for example, [H] [22]. Actually: 

Remark 1 ((Numerical) Taylor series method) Taylor series may be used 
as a tool to numerically solve the initial value problem associated with (TJJ). To 
this end the convergence of the Taylor series must be controlled. The most 
frequently used procedure is the stepwise (or piecewise) procedure described in 
section Wl.2.1\ 

In this respect of the numerical treatment of ODEs, and considering the high 
level of numerical development of the traditional Taylor method (see section 
I2.1.2.3p . one is forced to acknowledge that the current use of the DTM cannot 
be seen as a real technological break-through (despite the huge number of 
publications on the subject). One may easily verify that, the DTM (as it is 
currently used) brings no new capability into the numerical treatment of (at 
least ordinary) differential equations (see section [2]). Consequently, to possibly 
distinguish it from the traditional Taylor series method only remains the way 
the expansion coefficients are calculated, in apparent accordance with the 
above quotation of [T5j . 

Unfortunately, when looking at the "old papers" using the Taylor method, in 
that respect of calculating the terms of the Taylor expansion, one quickly real- 
izes that the basical tricks of the DTM (that transform the differential equation 
into an algebraic iteration equation, see section [2]) were already known and 
used (for a long time, indeed, see section I2.1.ip without having any need to 
change the name of the Taylor series method (e.g., see [T8lll9f21] ). Even some 
efficient procedures of the traditional approach have remained ignored by the 
DTM users (see section I2.2.2p . Several articles that present numerical soft- 
wares using the traditional Taylor series method with a view to automatically 
solve systems of ODEs are still currently published [23] -[33]. Based on an old 
and well established method (see, e.g. [32]), the authors of these articles have 
no need to mention the existence of the DTM, though they use exactly the 
same process! In fact, the reverse should be the rule: the DTM users should 
cite traditional studies, what they never do, probably because they do not 
know their existence. 

The status of the DTM is made very confused because that aspect of the DTM 
to be potentially a computer-specialized procedure is often ignored by several 
DTM users. Most authors using the DTM do not implement its numerical 
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ability, they provide the solutions to the problems studied under merely raw 
Taylor series. When such studies are directed towards systems of pure ODEs, 
the original contribution of the DTM is practically vacuous since it obviously 
coincides with extremely simple, to not say trivial, traditional studies (see 
section [2]) . 

This explains and justifies the recent biting comment [8] in which one may 
read: "the DTM is the familiar Taylor series method with a fancy name" and 
also: 

"There is no doubt that the DTM is merely the textbook Taylor series ap- 
proach that some authors have disguised as a fancy discrete transformation." 


The sharpness of the misunderstanding is well emphasized by the reply [H] to 
that comment: 

" Several revered authors and reputed journals [would] have nurtured this 'mis- 
conception ' over the years. It [would be] a challenge to such authors and jour- 
nals who have published several works on DTM." [5] 

It is time to clarify the status of the DTM. This is the tentative object of the 
present article. 

As one may have already foreseen, there are two aspects of the DTM, it may 
be seen either 

(A) as a " computer-specialized procedure" (as proposed in [1] for example), 
or 

(B) as a technique for determining the terms of the raw TayloiO] series of 
solutions of various problems (see [13] as an example among several oth- 
ers), 

or both, of course. 

Concerning interpretation A, it is worth recalling that, in general, the process 
of solving a differential equation using a Taylor (power) series requires two 
steps: 

(1) Obtening the solution as a formal power series (Taylor series, whatever 
the technique used to get it). 

(2) The eventual implementation of an analytic continuation process when 
the range of convergence of the direct sum of the series is too limited to 

5 More generally of power series including non-integer power series. 
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reproduce the solution sought in its expected domain of definition. This 
procedure may be used: 

(a) either, in a numerical routine to solve an initial value problem (or 
within a shooting method procedure) with fixed initial parameters 

(b) or in a semi-analytic approach of, e.g., a boundary value problem 
(the series is expressed in terms of the unknown -yet undetermined- 
initial parameters, see e.g., [33]). 

Now, it is true that, from its origin, the DTM implicitly involves the possi- 
bility of automatically implementing an analytic continuation via a piecewise 
process (although it is said in terms closer to the numerical version of item 
2a, as a multistage numerical method, see for example [I]). This possibility 
has been exploited by several authors (see [T2lll4|ll6|[TT| and [35] [50]) . In that 
case, the method is truly different from the raw Taylor series method which, 
stricto sensu, does not involve any consideration of analytic continuation (see 
for example [M]). One may then rightly consider the method as being "an 
extended Taylor series method" [H]. But, as explicitly shown in section [2j in 
the circumstances, the DTM merely uses the "analytic continuation" process 
of the Taylor series, as already Euler did it (see jH]) and which is described in 
textbooks under the name of "continuous analytic continuation" |12J chap. 9] 
or of "Taylor algorithm" |33l p. 267]. This is a technique included for a long 
time and still currently used in numerical softwares for solving ODEs by the 
"traditional" Taylor series method (see, [H] [33] and [4~4"f4"5] ). 

There is a more serious problem with the way that the DTM is viewed by some 
authors. This problem is linked with the interpretation B and the analytic 
continuation presently in question. 

Many studies of ODEs using the DTM have been applied according to inter- 
pretation B to so simple problems that the solutions sought were analytical 
(at least in the range of supposed interest of the independent variable, it is, 
in particular, the case of the study criticized in [8]). In such circumstances, 
where no analytic continuation was needed or performed, one has employed 
the expression DTM (the name of the technique used to get the raw Taylor 
series) as if it was a new (formal) method in place of the Taylor series method 
itself. This unfortunate view of the DTM is, by now, so common that when 
an analytic continuation process of the series has been realized by means of 
Pade approximants, one has changed the name of the DTM into the modified 
differential transformation method (MDTM or also the "DTM-Pade method" ) 
|46 | |51j . Such names are not justified because the DTM enables, in princi- 
ple^, the automatic implementation of Pade approximants (see section I2.2.3p . 



6 It is possible (see footnote [2]) that this eventuality has been forseen in the original 
presentation of the method, but -to my knowledge- this has never been used by the 
defenders of the DTM in the current literature, hence the new name of MDTM. 
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Anyway, in the circumstances, the (formal) method used is merely the Taylor 
series analytically continued by Pade approximants (see, for example, [31]) a 
method which has been used in numerical routines for a long time |52j . 

Even more serious is the recent and repeated proposal of introducing, as a 
novelty, a piecewise procedure within the DTM and of naming this "new 
method" again MDTM (with "M" for modified or multistage or multistep) 
[53] [58]. The authors of such proposals (and the referees who have allowed 
their publications) were manifestly unaware of the fact that the procedure was 
already clearly foreseen in the original DTM (see, for example, PQ) and already 
implemented by several DTM users [T^|14|16|IT7] . [35] EO] , not to mention the 
previous utilizations (some are cited above) by the traditional-Taylor-method 
users. 

Extensions of interpretation B to problems which are potentially more tricky 
than a simple system of ODEs are seemingly more interesting (despite the 
standing formal character of this interpretation). As a matter of fact, the 
DTM has been used as a technique for determining the terms of the raw 
poweiQ series of solutions of various problems. Its basic ingredients, though 
not new in the current use, have been presented in such a manner that it has 
appeared very attractive (see section [272]) . Rid of the original constraining goal 
of performing finite precision calculations, its implementation, sufficiently easy 
and clearly presented, has been formally extended to several different kinds 
of equations. Thus, equations of the following types have been considered 
using the DTM: integro-differential [591160] . partial-differential [39], [6l]-[70]. 
differential-difference [51~|71f72] . differential-algebraic [?3]r74"] . fractional order 
differential [58], [75]-[90], fuzzy differential [91], q-partial differential [92], and 
so on. Of course, some of those subjects had already been similarly (i.e. an- 
alytically) addressed by Taylor-method users, but the lack of a systematic 
formalization of the traditional approach could have limited their impact (see 
[93] [103] . for purely analytical approaches). However, if one considers studies 
with finite precision goal then, again, the traditional method has reached a 
higher level of development such as in the treatment of the differential alge- 
braic case -a particular system of ODEs already well treated [2"3"f3"2"] . [104] - 
[107] - or of integral and integro-differential equations (see, e.g., [96|108j ). I will 
not discuss further the use of the DTM in such problems. Nevertheless, it seems 
that the DTM could deserve its name (as a technique) when it extends the 
Taylor method to new kinds of expansions (different from a Taylor expansion, 
see e.g., [109] ) and to unusual derivation processes such as in the q-calculus. 
Nevertheless, in such studies, considering their formal character limited to the 
obtaining the algebraic iteration equations for the coefficients of the series ex- 
pansion, it should be clearly stated that the method used (to solve the primary 
equation) remains merely the Taylor (or power) series method. 



Including non- integer power series. 
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Moreover, with interpretation B, the series is calculated only once (and most 
often it is limited to low orders). Consequently the recourse to an iterative 
procedure to determine the raw power series is even not necessary. Thus, in 
this meaning, the recourse to DTM may appear entirely artificial (a "fancy" 
name). As already said, that version of the DTM is so commonly used that 
it is by now often presented as "the" DTM. That interpretation, surely, has 
induced the polemic and the confusions mentioned above. This does however 
not give more credit to interpretation A, as shown in what follows. 

The following section is limited to the problem of solving ODEs (numeri- 
cally or not). It is explicitly shown that the traditional Taylor method should 
take precedence over the DTM. It is also indicated why it is not acceptable 
to modify the name "DTM" when the resulting Taylor series is analytically 
continued. 



2 The Taylor series method vs. the differential transform method 

In this section, a comparison is made between the most elaborated current 
formulations of the Taylor series method and of the DTM when they are 
applied to solve an initial (or boundary) value problem of ODEs (for recent 
reviews on the use of the Taylor method see, e.g., [21] [33] and [HOj and for 
the most elaborated presentations of the DTM see, e.g., [TOl2|14|40|lll] ). 
Only the one dimensional case of the first order ODE ([T]) is considered here, 
the generalization to several variables and to systems of ODEs of higher orders 
is straightforward. 

2. 1 The Taylor series method 

The recourse to "infinite series" to solve differential equations is an old method 
which, as mentioned in jHJ p. 4], may be attributed to Newton in [1 12j . Since 
then, the method has been rediscovered several times (as well illustrated, 
notably, in |113j ). By now, it is named the Taylor series method. Several 
current papers may be found that present this method in details (for example: 
[32|33|110j . see also jJH pp. 47-49]). The following considerations are thus far 
from being original. 

For the sake of considering practical solutions of flTJ), and with a view to 
describe the "traditional" Taylor method, definition [TJ given p. [3] must be 
slightly modified. 

With the Taylor series method, the generic solution of ([T]) is expressed locally 
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truncated Taylor expansion about a time U: 



x(t) = Ek=o (t-t i ) k X k + 0[((t-t i 
k = 1,2,- •• ,N. 



N+l 



(3) 



Xk — — 



l_ d k x 

k\ HF 



t=U 



In the following it is assumed that t > t i} and a, is supposed to be known 
either because ti is the initial time (say to) & t which the initial condition is 
provided or because it is an intermediate time value (ti = t + ih) reached by 
a stepwise analytic continuation procedure (see section [2.1. 2. ip . Although not 
indicated here and in the following (except where noted), the coefficients Xk 
depend on ti. The dependence is implicit via the initial value oti, but it may 
also be explicit if the ODE considered is non-autonomous. 

As said in the introduction, the method involves two steps: first the obtaining 
the series for some (higfu 8 1) value of N, and second a process of analytic 
continuation of the series. Stricto sensu, the second step is not a part of the 
formal Taylor series method. But, because it is by now difficult to distinguish 
between studies based on Taylor series which are formal (getting only the raw 
series) or numerical (using, in addition, analytic continuations), I temporarily 
include the two steps in the traditional Taylor series method (awaiting to show 
below why that is justified), reserving the name raw or formal Taylor series 
method exclusively to the first step (in accordance with definition [T] p. [3]). 



2.1.1 Getting the raw Taylor series 



The coefficients Xk must be determined such that the ODE flTJ be formally sat- 



one condition^ (x (ti 



isfied order by order up to O ((t — U 
quired so that the resulting system of N equations for the Xk (k 
has a definite solution. 



= cui) is re- 
0, N) 



Establishing and solving that system of equations is not the most efficient 
way of determining the coefficients Xk (in particular this would imply the 
explicit calculation of "the necessary derivatives of f(x,t)"). In fact, it is 
known for a long time by the (thus traditionnal) Taylor-series-method users 
that the Xk may be determined more efficiently by iteration provided / (x, t) 
be simple enough (see, e.g. |114|I115|I116] . see also [UJ pp. 47-49]). Instead 



8 Actually, it is rarely necessary to consider high values of N, as is shown in the 
following. 

9 In general the order is N — p, with p the order of the ODE, and p conditions are 
required to unambiguously determine the X^. 



9 



of N equations, only one equation (the iteration equation) has, then, to be 
effectively considered. 

Assuming, as a simple illustration, that f (x,t) = Xx(t) with A a constant, 
it is easy to show that the generic Taylor coefficients X k satisfy the following 
iteration equation: 

X = OLi , 

(k + 1) X k+ i = \X k , k = 0, 1, ■ ■ • , oo . 




This is a difference equation of same order (first) as the ODE under study. 
With the initial condition, the coefficients X k are obtained easily by iteration 
of (jlj) up to high orders, with e.g.: 

A 2 A^ 
X = ai, X 1 = a i \, X 2 = at— , X N = at— , 



The terms of the Taylor series of an exponential is recognized giving the exact 
solution (which satisfies the right initial condition): 

x(t) = ai e x{t - u) . 

Although that example is trivial, it indicates the way towards a more elabo- 
rated procedure when / (x, t) is more complicated. The above example may 
be presented as follows. One constructs a priori a table of correspondence: 

x (t) -> X k , 

-^(k+l)X k+1 , 
Xx(t) -> \X k , 

so that the ODE of interest is automatically transformed into a difference 
equation: 

^ = Xx(t)^(k + l)X k+1 = XX k . 

Generalizing to the k th Taylor coefficient of / (x,t) (denoted below by F k ) 
then it comes: 

{k + 1) X k+1 = F k . (5) 

To construct the F k s, the Taylor-series- method users have considered supple- 
mentary obvious (or already known) correspondences such as (see [1171 p. 525] 
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but also [Tig] ): 



(k + n)--- (k + l)X k+n , 
x(t) + y(t) ^X k + Y k , 
x{t)y(t) -+I$ =0 X i Y k - i , 



(6) 



* (*) = if i (** - Et-o 1 ^n-i) , 

«(t) = [x(t)]^E^i(^i-l)^*-i > 



fc > 1, Z = X, 



One has then considered more complicated forms of / (x, t) such as exponential 
(similar to [x(t)] ), logarithm, trigonometric functions etc... |18|ITU] (see also 
[JU pp. 48, 49]) so as to find similar simple formulas for their images. The 
aim was, notably, to include them into numerical softwares |119j with a view 
to solve automatically various kinds of systems of ODEs [T8|19|ll20j . 




As noted in [22], such recurrence schemes for the terms of a Taylor series in 

solving a differential equation was used as early as 1946 by J. CP. Miller |114j . 
Some authors (see, e.g. [121] ) attribute the first use of such recurrences to 
the computation of derivatives by J.R. Airy in 1932. Actually, the iterative 
procedure was known already at the time of Newton (see, [1 12] Prob. II, Sol. 
Case II, Ex. I, p. 33] and also, e.g., |122[ p. 116]), and the last line of (jSJ), often 
attributed to J. CP. Miller |117[ p. 507 of second ed.] and [123] , had been 
established by L. Euler in 1748 [124] . cf [1171 P- 526]. Such a formula is the 
result of the application of the Leibnitz rule for the derivative of a product 
of functions. This rule, and others, is at the basis of an efficient procedure 
for numerically calculating high order derivatives named automatic differenti- 
ation (for extensive bibliographies, see [125|I126] ). which is used in turn to get 
efficient iterative formulas for complicated forms of f(x,t). This technique is 
used for a while in numerical routines [23], [22] -[HT], [331121] . [127] - [T3U] and 
in regular studies (e.g., [IS]). More recently, similar considerations have led to 
a proposal [131] for efficiently implementing the Picard iteration method. 

It is worth indicating that such recurrences (and others) have been established 
also for functions of two variables (see e.g., [20] for analytical expressions) 
and that the automatic differentiation process in many variables is also well 
developed for numerical applications (see, e.g., |132|I133|I134] ). 

In the following, the name "Taylor transformation" will designate the above 
described automatic procedure for obtaining the iteration equation (the image) 
of a given ODE (by the way, it was the primary name of the DTM p] a name 
rightly re- used in the first papers of the 90's [T0|11II12] ). 




It is important to keep in mind that, up to now, the above rules of calculation 
of the terms of the Taylor series is purely formal. In particular it does not an- 
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ticipate the nature of computation (theoretical exact or finite precision) which 
is to be done. In the following subsections, the considerations are less formal 
and are progressively more and more oriented towards numerical applications. 

2.1.2 Analytic continuations 

Sometimes the generic Taylor series is convergent in the domain of definition of 
interest (such as in the preceding trivial case of the exponential). The (formal) 
method then provides an exact solution (when N = oo and if the generic term 
is identified) or, at least, a convergent form towards the exact solution (for 
iV finite). Those simple cases, often linked to linear ODEs, are not of great 
interest. 

More interesting are the nonlinear ODEs, the generic Taylor series of the 
solutions of which have finite radius of convergence due to the presence of 
singularities in the complex plane of the independent variable t. The efficiency 
of the Taylor expansion thus is limited to a finite domain (which may be small) 
about the expansion point tj. To enlarge this domain, the recourse to analytic 
continuation procedures is often required. (Notice however that in practice, 
one may need to use analytic continuation even for series with an infinite 
radius of convergence, if only to reduce the length of the series.) 

There are several kinds of analytic continuations. Let us consider only the two 
most often used in conjunction with the Taylor series method. 

2.1.2.1 Stepwise procedure The stepwise procedure, also named the 
"continuous analytic continuation" [121 chap. 9] or the "Taylor algorithm" [13| 
p. 267], is particularly interesting because, beyond its character of analytic 
continuation procedure, it is structurally well adapted to numerical routines. 

Let the domain of definition of the solution sought, say t G [to, to + H], be 
cut innGN pieces: H = nh. Then the times of reference read tj = to + ih, 
% — 0, 1, . . . ,n. If h is chosen sufficiently small (at each step, smaller than the 
distance to the closest complex singularity of the solution), then the Taylor 
series converges in the domain [ti,t i+ i], and its sum provides the solution in 
this domain with a sensible accuracy even if N is not very large. If this is 
true for each value ti then one gets, step by step, an approximate solution in 
the whole domain [0, H], each initial value CKj being (approximately) provided 
by the sum of the Taylor series at the second boundary of the preceding sub- 
domain, namely 

JV 

a i = , Eh k (X k ) ii _ 1) + 0(h lf+1 ) , (7) 

fc=0 
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in which I have indicated the dependence of X}. on the step % — 1. 



This stepwise procedure has been (and is currently) successfully used in nu- 
merical routines by Taylor-series-method users to solve initial value problems 
(for reviews see, e. g. 1241132(33] ). Its first use in conjunction with power series 
goes back to Euler [UJ. For that reason, it is justified to consider it as being 
an integral part of the traditional Taylor series method as the name " Taylor 
algorithm" (131 p. 267] suggests it. 

The procedure has also been used to solve boundary value problems [T6"f23|ll20|n 
In that case, the unknown «o of the initial boundary is transmitted, step by 
step via the parameters oti, up to the second boundary t + H where the de- 
sired condition is imposed. The value of the unknown is chosen among the 
acceptable solutions of the resulting equation. The method is limited by two 
aspects: 

(1) The dependence of oti on the unknown becomes rapidly very compli- 
cated. (However, see |136|ll37j in which supplementary Taylor expansions 
in powers of the initial unknown are used as generally done when using 
the Taylor method in a sensitivity analysis, see e.g., [133|I134] .) 

(2) When H — > oo the method no longer works because, when t — > oo the 
truncated Taylor series goes to ±oo according to the sign of the last term 

In the last case, where the domain of definition of the problem is infinite, it is 
better to consider other kinds of analytic continuation such as, for example, 
rational fractions (Pade approximants) . 

2.1.2.2 Pade approximants The use of Pade approximants is well known 
[138|ll39|ll40j . and will not be described here. It suffices to recall that it is a 
procedure that transforms a truncated series (a polynomial) into a ratio of 
polynomials. It is a so drastic modification of the Taylor series of the generic 
solution that it cannot be easily seen as the "inverse" of an "image" resulting 
from some "modified" Taylor transformation. However, it is shown in sec- 
tion 12.2.31 that, similarly to the step- wise procedure, the recourse to Pade- 
approximant seems to have been also foreseen in the DTM but, seemingly, 
has never been implemented in the current literature. 

2.1.2.3 High level numerical development There is no doubt that, as 
a numerical tool, the traditional Taylor series method has reached a high level 
of developement. I have already mentioned several numerical routines that 

10 Euler 's motivation was not to perform an analytic continuation but, precisely 
in opposition to the statement of |15j quoted in the introduction, to avoid the 
calculation of high-order derivatives of / (x,t). 
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have been developed with a view to automatically solve the initial (and bound- 
ary) value problems of systems of ODEs [23] [33]. I must mention also the exis- 
tence of adjustments of the Taylor method to treat stiff systems [23].[25]-|52], 
|141j — and the possibility that this method offers to get high-precision so- 
lution for ODE |146j . In the same vein, it is worth mentioning the use of Taylor 
series methods in interval (or validated) solutions, that provide a garantee of 
existence of (and bounds on) the solution in a given interval |129|ll30j . |147j - 
|151] (for a bibliography on the subject of validated solutions see, e.g., the 
paper by Nedialkov et al in |129|ll3Uj ). As an extension to those latter stud- 
ies, it is fair to mention also the developement of new methods for obtaining 
garanteed bounds on the expansion |136j . |152j — p~&4~] . for a review see Makino 
and Berz in |155j . 

In this respect of the quality of the numerical treatment of ODEs the differ- 
ence is glaring with studies using the DTM in which almost never are raised 
questions like, e.g., selecting the stepsize (or the order of development) or dis- 
cussing the computational complexity of the system considered. In particular, 
despite the original Pukhov's aim, I have not seen any paper that proposes a 
numerical routine with a view to automatically solve systems of ODEs nor a 
possible adaptation of the DTM to treat stiff systems. 



2.2 The differential transform method 

In this section the most elaborated versio of the DTM and its usual version 
are presented and compared with the Taylor series method presented in the 
preceding sections. 



2.2.1 General procedure 

The DTM [I] -[5] is a formalized modified version of the Taylor transformation 
described in section 12.11 The basic difference is that instead of Xf~ defined in 
©, one considers the image X k of x (t) defined as [TTTIT2TT4T4WTT] : 

y d" (g(t)x(t)) 

in which the auxiliary function q (£) and the infinite set are chosen a priori. 
Mfc is called the weighting factor and q (t) is named the kernel corresponding 
to x (t) (see, e.g. jlllj). 



Accessible to me, see footnote [2j 



t=t, 
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One recovers x(t) through the following series (the inverse transform): 

* E^-^#- (9) 



2.2.2 Usual procedure 

To my knowledge (see footnote [2]) , the freedom in the choice of q (t) has never 
been exploited in the current literature, where it is usually set equal to one 
(for possible uses of q(t), see section 12.2. 3. ip . Moreover, the choice of M k is 
currently limited exclusively to the two following cases (except in jlD] where 
an adaptive grid size mechanism is employed): 

• M k = |f, then X k = X k (since q(t) = 1), the DTM coincides with the 
Taylor transformation and the method (to solve a differential equation) 
to the formal Taylor method (without any analytic continuation -i.e., the 
interpretation B mentioned in the introduction). 

• M k = as proposed in [IJ3] . where h has the same meaning as in section 
12. 1.2. H and is effectively used in the same way (see, e.g. [3]). One has X k = 
h k X k , and eq. (J7J) is slightly simplified: 

N 



fc=0 



The drawback is that one must account for the presence of h in the defini- 
tions of the images such as, e.g. 

x it) -> X k , 

dx(t) . (fc+1) y- 
~dt * ~h~ A k+l ■ 

It is easy to verify that one obtains the same result by merely rescaling the 
time variable t —¥ ht within the traditional Taylor approach. Although, in 



practice, there may be advantages one way or the otheii 12 1, this is too little 



to justify a drastic change of designation of the method used. 

For the remaining aspects of the DTM, that is to say the way that / (x,t) is 
treated, there is no difference from the Taylor transformation described in sec- 
tion [27TTTJ Only simple forms of / (x, t) are considered, and each case is treated 
separately. Finally a table of correspondence similar to (jfJJ) is constructed. 



12 Indeed, this form has been included in some routines of Taylor-series-method 
users to reduce underflow and overflow in practical computation by choosing h of 
the order of the radius of convergence of the series [22] , 
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It is important to mention that, concerning the construction of the table of 
correspondence, DTM users have recently been reminded of the existence of 
the exponentiation recurrence [last line of (JO])] with a view to improve the 
efficiency of their method [156] . The same kind of remark had already been 
made some time before in |157|I158] where the recourse to the Leibnitz rule 
of derivation is used to get the image of the exponential or logarithm func- 



tions through "new algorithms! 13 1 which are more efficient than those of the 
standard DTM. These proposals for improving the DTM clearly contradict 
the claim quoted in the introduction that the traditional Taylor series method 
would be "computationally intensive as the order becomes large" compared to 
the DTM since, for a long time, those recurrences were known (see, e.g. |118] 
and section l2.1.ip and already included in routines to solve ODEs by Taylor se- 
ries (see, e.g. [20|22] ). These late propositions for improving the DTM, and the 
succes of the DTM itself, show enough that, as the traditional Taylor method, 
the automatic differentiation procedure |125j - [T3"U] is largely unrecognized in 
the current litterature, this is a pity. 



2.2.3 The modified differential transform methods 

2.2.3.1 Modification by Pade approximants When the domain of 
definition of the solution sought is infinite (H = oo), some authors have pro- 
posed substituting for the stepwise procedure foreseen in the DTM by using 
Pade approximants [16] [51]. To this end one first chooses = whereas 
q (t) is still set equal to one (i.e. the raw Taylor series is obtained first). One 
then transforms the Taylor series (the inverse image) into a rational frac- 
tion (according to the rule of Pade calculus [138] - [T3D] ). Now, because the 
recourse to Pade approximants (eventually after a Laplace transform of the 
series [15911160] ) has been implemented at the level of the inverse of the image 
of x (t), then this is no longer a differential transformation (modified or not). 

In fact, eq. (Q shows that the DTM enables the automatic implementation 
of the Pade-approximant technique. Indeed, it suffices to decide that q (t) is 
a polynomial in (t — ti) of order N 2 (with coefficients Qk, k = 1, • • • , N 2 to 
be determined and Qo = 1) and to limit the X^-series to be of order N\ with 
Ni + N 2 = N, to get a balanced system of equations for the Qk and Xk 
expressed directly in terms of the "image" of / (x, t) and of the initial value 
aj. To be concrete, the system of equations for the Qk's and X k 's is: 

Y^ O Q i X k - i = ) fc = JVi + l,... ,N X + N 2 , 

Xk = Si=So Qi Xk-i, k = 0, • • • , Ni , 
in which Q Q = 1 and the X^s are implicitly known from §5$) in terms of the 
13 Already known by the users of the traditional- Taylor-series method. 
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initial value ctj. 



Notice that, because the first line of ( TTTj) is a system of N 2 linear equations 
which cannot be solved iteratively, it is necessary to determine the X^s explic- 
itly Consequently, the process is not more efficient than applying the Pade- 
technique on the Taylor series and, in any way, the method used remains 
merely the Taylor series analytically continued by Pade approximants. 



An other possible usd 14 1 of q (t), in the case of an infinite value of the range 
H, could be as follows. Sometimes one knows that, when t — > oo, the solution 
sought behaves as f with v e Q, and this behavior is hardly reproduced by 
a Taylor series even continued by a Pade approximant. To circumvent this 
difficulty, one may choose (assuming ti = 0): 

q(t)=(l+ t -) U (12) 

in which r is a parameter to be adjusted in order that the radius of convergence 
|r| of the series expansion of f fl2|) be larger than that (a priori unknown) of 
the Taylor series of x (t). With such a choice for q (t), the behavior at infinity 
of the transformed function is a simple constant reproducible by a diagonal 
Pade approximant. 



2.2.3.2 Modification by the stepwise procedure I have already ex- 
pressed, in the introduction, what must be thought of the proposal of improv- 
ing the DTM by introducing a stepwise analytic continuation of the Taylor 
series, when this procedure was already clearly foreseen in the DTM and used 
several times as recalled in section 12.2.21 



2.2.4 Formalization 

With a bit of emphasis, the interpretation B of the DTM has been presented 
on the same footing as the Laplace, Fourier, and Mellin transforms where 
the integration procedure used to construct the image of a given function is 
replaced by a derivation procedure [5] (see also |161j ). The inverse transforma- 
tion gives a power series (when q (t) — 1 and — it is merely the Taylor 
series of the function). By similitude with the Laplace etc. transformations, 
a list of properties of the DTM with respect to the addition, multiplication, 
etc. has been given, as well as the images of current functions. The result- 
ing transformation of, say, an ODE into an algebraic equation easily solvable, 
has been presented as being formally equivalent to that of the Laplace (or 

14 Notice also that could be used to implement a Borel transform of the Taylor 
series of x (t) . 
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Fourier) transform with the additional advantage that the DTM applies also 
to nonlinear ODEs. 

Of course, the depth of the considerations is lesser than in the traditional 
Taylor series approach. But this latter has been almost exclusively developed 
as a numerical tool. This ability of the Taylor series method being not very 
well publicized, the formal presentation of the DTM has appeared (wrongly) 
as a new and elegant way of constructing a table of correspondence like fl6]). 

As emphasized enough in the present article, given the high level of develope- 
ment of the traditional Taylor series method (as a numerical tool), the DTM 
contributes nothing new in the way systems of ODEs are solved. Nevertheless, 
the DTM has allowed an easy generalization of the Taylor method to various 
derivation procedures. For example, fractional differential equations have been 
considered using the DTM extended to the fractional derivative procedure [58] , 
[ZH]-[9U| via a modified version of the Taylor series |109j . In similar situations 
-such are also fuzzy differential [91], q-differential [92], etc.- the DTM might 
deserve its name (although traditional Taylor series treatments also exist that 
deal with fractional derivative |103] or fuzzy equations |101|102j ). 

Less evident is the original contribution of the DTM to the study of partial 
differential, integro-differential, integral, etc. equations where the usual deriva- 
tive procedure is in action (e.g., integral and integro-differential equations have 
been studied by Taylor series methods in [93]-[98], |108j ). 

In general, it would be fair to not forget the limit of the formal version of the 
DTM which is not a method to solve the initial problem but only a conve- 
nient way of obtaining the iteration equation for the power series coefficients 
(provided the rules of the automatic derivation be well accounted for). For the 
sake of clarity, in such situations, a clear reference to the traditional Taylor 
(or power) method should be mandatory. 



3 Conclusion 

A misunderstanding on what the "traditional" Taylor series method is has 
grown out since the 90's of a desire to promote an attractive method named 
the differential transformation method. Actually, when the DTM was being 
born, the Taylor series method had developed for quite a while the treatment 
of ODEs exactly on the same basis and, concerning this treatment, absolutely 
nothing new has originated from the new method. Even some efficient recur- 
rence formulas (known for a while) have been unknown to the DTM users. 

The misunderstanding has increased over time because the DTM (and also 
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the Taylor series method, of course) has two possible uses: as an analytic tool 
or as a numerical tool. Although it was created with this latter use in mind, 
the DTM has most often been used as an analytic tool, so that some authors 
have even forgotten the initial aim and proposed to reinvent the method. 

Beyond the treatment of ODEs with the DTM -in which case one should at 
least clearly refer to the numerous studies done with the Taylor series method- 
it seems that the major contribution of the DTM is in the easy generalization 
of the Taylor method (either analytical or numerical) to problems involving 
unusual derivative procedures such as fractional, fuzzy or q-derivative. 
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